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ABSTRACT 


One of the primary goals of the NASA Sensor Intercomparison and Merger 
for Biological and Interdisciplinary Oceanic Studies (SIMBIOS) project is to de- 
velop methods for meaningful comparison and possible merging of data prod- 
ucts from multiple ocean color missions. The Modular Optoelectronic Scanner 
(MOS) is a German instrument that was launched in the spring of 1996 on the In- 
dian IRS-P3 satellite. With the successful launch of NASA's Sea-viewing Wide 
Field-of-view Sensor (SeaWiFS) in the summer of 1997, there are now two ocean 
color missions in concurrent operation and there is interest within the scientific 
community to compare data from these two sensors. In this paper, we describe 
our efforts to retrieve ocean optical properties from both SeaWiFS and MOS us- 
ing consistent methods. We first briefly review the atmospheric correction, 
which removes more than 90% of the observed radiances in the visible, and then 
describe how the atmospheric correction algorithm used for the SeaWiFS data 
can be modified for application to other ocean color sensors. Next, since the re- 
trieved water-leaving radiances in the visible between MOS and SeaWiFS are 
significantly different, we developed a vicarious intercalibration method to re- 
calibrate the MOS spectral bands based on the optical properties of the ocean and 
atmosphere derived from the coincident SeaWiFS measurements. We present 
and discuss the MOS retrieved ocean optical properties before and after the vi- 
carious calibration, and demonstrate the efficacy of this approach. We show that 
it is possible and efficient to vicariously intercalibrate sensors between one and 
another. 



1. Introduction 


The German Modular Optoelectronic Scanner (MOS) [1] is an imaging push- 
broom CCD spectrometer that was launched in the spring of 1996 on the Indian 
IRS-P3 satellite. The IRS-P3 satellite is in a sun-synchronous polar orbit with a 
mean altitude of approximately 817 km and a local 10:30 AM descending node 
crossing time, and completes one orbit about every 101 minutes. MOS has 13 
spectral channels ranging from 400-1010 nm with bandwidth of about 10 nm. 
The MOS image has a ground resolution of ~0.5 km in nadir viewing with scan 
swath of -200 km. MOS is a technology demonstrator instrument with limited 
geographic coverage capabilities. Its scientific applications are mainly in ocean 
color and atmospheric aerosol studies. With the successful launch of NASA's 
Sea-viewing Wide Field-of-view Sensor (SeaWiFS) [2] on August 1 of 1997, there 
are now two ocean color missions in concurrent operation. Therefore, we have 
an unprecedented opportunity to compare ocean color data from two sensors in 
simultaneous operation on two different satellite platforms. SeaWiFS has 8 
spectral bands covering from 412-865 nm with bandwidth of 20-40 nm. SeaWiFS 
is in a sun-synchronous polar orbit of 705 km altitude with a local noon crossing 
time and completes one orbit every 99 minutes. The SeaWiFS image has a 
ground resolution -1 km with a scan swath of -2800 km, enabling it to provide 
daily global coverage. As shown in Table 1, MOS has a slightly different spectral 
band characterization in comparison with SeaWiFS. It should be noted that only 
the MOS bands which are close to the SeaWiFS spectral channels are listed in Ta- 
ble 1. Also, the MOS band number is named corresponding to the SeaWiFS's as 
shown in Table 1 and used in this paper. The object of this paper is to develop a 
vicarious calibration approach in which the MOS spectral bands can be re- 
calibrated from the SeaWiFS measurements, thereby allowing remotely retrieved 
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ocean color results from the two sensors to be meaningfully compared. The vi- 
carious calibration method is also applicable for re-calibrating the satellite sensor 
with in situ ocean and atmospheric optical property measurements. 

2. Applying the SeaWiFS Atmospheric Correction Algorithm to MOS 

In this section we first briefly review the SeaWiFS atmospheric correction al- 
gorithm and its implementation into the SeaWiFS data processing system. Next, 
we present the modifications required to implement this algorithm for alternate 
ocean color sensors, e.g., the MOS. We then test the accuracy of the correction 
algorithm at the MOS spectral bands for various cases using the current SeaWiFS 
aerosol lookup tables. Finally, we compare retrieved results from MOS and 
SeaWiFS measurements using a consistent atmospheric correction algorithm for 
scenes acquired at various locations and different times. 

2.1 The SeaWiFS atmospheric correction algorithm 

We begin with definition of the reflectance p = kL//j. q F 0 , where L is the radi- 
ance in a given solar and viewing geometry, F 0 is the extraterrestrial solar irradi- 
ance, and p 0 is the cosine of the solar zenith angle. The total reflectance meas- 
ured at the top of the ocean-atmosphere system can be written as: 

PAX) = pAX) +P.(A)+ p„(X) + i{X)p wc {X) + i(X)p.(X), (1) 

where p r (X) is the reflectance resulting from multiple scattering by air molecules, 
p a (X) is the reflectance resulting from multiple scattering by aerosols, p ra (X) is the 
multiple interaction term between molecules and aerosols [3], p wc (X) is the re- 
flectance at the sea surface that arises from sunlight and skylight reflecting from 
whitecaps on the surface [4], and p w (X ) is the water-leaving reflectance, which is 
the desired quantity in ocean color remote sensing. The t(X) is the atmospheric 
diffuse transmittance [5, 6] that accounts for the effects of propagating water- 
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leaving and whitecap reflectances from the sea surface to the top of the atmos- 
phere (TO A). In the above equation, the surface sun glitter term has been ig- 
nored. Observations which have significant sun glitter contamination can not be 
accurately corrected, and must be removed. To relate the derived water-leaving 
reflectance to the ocean inherent optical properties (IOP), the atmospheric effects 
on the water-leaving reflectance p w (X) must be removed. The normalized water- 
leaving reflectance, [p vv (A)] yv , can be defined from Gordon and Clark [7], 

(P.a)]» = P»(A)//(A,0 o ), (2) 

where f(A,0 o ) is the atmospheric diffuse transmittance in the solar direction with 
the solar zenith angle of 6 0 . The value of two-band ratio of [p w {X)] N in the visi- 
ble can then be used to infer the ocean near-surface optical properties [8-10]. 
Note that in comparing the retrieval results from two different sensors which 
usually have slightly different spectral band characterizations, the normalized 
water-leaving reflectances provide a more meaningful comparison than the radi- 
ance values. The radiance value is a function of the solar irradiance, which will 
vary with a sensor band's spectral response. 

Since more than 90% of the signal in visible measured at satellite altitude is 
contributed by the atmosphere (the first three terms in Eq. (1)), accurately re- 
moving the atmospheric effects is crucial to the success of any ocean color remote 
sensing experiment. The Gordon and Wang atmospheric correction algorithm 
[11] uses the SeaWiFS two near-infrared (NIR) bands centered at 765 and 865 nm 
to estimate the atmospheric effects and extrapolate these into the visible. Unlike 
Rayleigh scattering which can be computed accurately, the aerosol scattering is 
highly variable, and the effects of the p a (X) + p ra (X) in Eq. (1) on the imagery 
cannot be predicted a priori. The water- leaving reflectance p w ( X) at the two NIR 
bands, however, is usually negligible because of strong water absorption. There- 
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fore, the radiances measured at these two NIR bands are essentially the con- 
tributions from the atmosphere. For the SeaWiFS two NIR channels, Eq. (1) can 
be written as 

p,(A) - p r ( A) - t{X)p wc (X) = p a (A) + p rfl (A). (3) 

Therefore, the effects of aerosols and Rayleigh-aerosol interactions, p a (A) + p ro (A), 
in the imagery can be estimated at the two NIR bands from the sensor-measured 
radiances, the computed Rayleigh scattering reflectances, and the estimated 
whitecap contributions [4]. This quantity is then extrapolated and removed in 
the visible. The extrapolation was achieved through a process of aerosol model 
selection from evaluation of the atmospheric correction parameters, e(i,j), de- 
fined as [11-13] 

£(i’j) = Pas(i)/PasU)' (4) 

where p^i) is the single scattering aerosol reflectance at a wavelength A, . The 
Xj is usually at the longer NIR band, i.e., 865 nm. The value of e(i,j) character- 
izes the spectral variation of aerosol optical properties which include the aerosol 
optical thickness, single scattering albedo, and the aerosol scattering phase func- 
tion. It therefore can be used to infer the aerosol models. 

The implementation of the Gordon and Wang algorithm into the SeaWiFS 
data processing system was achieved through the use of lookup tables based on a 
large number (-25,000) of radiative transfer simulations that use the 12 aerosol 
models developed by Shettle and Fenn [14]. The main lookup tables contain in- 
formation of the p a (A) + p ra (A) values for various aerosol optical and microphysi- 
cal properties (12 aerosol models with various aerosol optical thicknesses) and 
solar and viewing geometries at the 8 SeaWiFS spectral bands. Generating the 
aerosol lookup tables involves a large number of radiative transfer simulations 
and requires substantial computer resources. 
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2.2 Atmospheric corrections for MOS 

Application of the SeaWiFS atmospheric correction algorithm to MOS would 
be difficult if it was necessary to regenerate the aerosol lookup tables for the 
MOS spectral bands. In a recent paper, Wang [15] discussed the effects of spec- 
tral band variation on the SeaWiFS atmospheric correction algorithm and out- 
lined simple procedures necessary to implement the algorithm for other ocean 
color sensors. In summary, to apply the SeaWiFS atmospheric correction algo- 
rithm to MOS we need to, according to the MOS spectral response functions: (i) 
re-compute the extraterrestrial solar irradiances and ozone absorption coeffi- 
cients, (ii) re-generate the Rayleigh scattering radiance tables at the sensor's 
spectral bands, and (iii) modify the atmospheric diffuse transmittance computa- 
tions. Of these steps, procedure (ii) is the most important. 

We have implemented the SeaWiFS atmospheric correction algorithm for 
MOS and tested algorithm performance for the MOS spectral bands with simu- 
lations. Following Gordon and Wang [11], we have applied the correction algo- 
rithm to a series of simulations carried out using the Maritime aerosol model 
with the relative humidity (RH) of 80% (M80 refers to the Maritime aerosol with 
RH = 80%), i.e., p t (X) was simulated with M80 aerosol model at the MOS spectral 
bands assuming that pjX) = 0. The SeaWiFS aerosol lookup tables, p a (X) + p ra (j l), 
were used for all computations. The error in the retrieved water-leaving reflec- 
tance, Ap(X) = t(X)Ap w (X), was computed. Fig. 1 provides results of algorithm per- 
formance for the MOS spectral bands at different solar and viewing geometries 
for the M80 aerosol model with aerosol optical thickness of 0.2 at 865 nm. For 
reference, a 5% error in water-leaving radiance at 443 nm, which is the SeaWiFS 
goal, corresponds to Ap -0.001-0. 002. Fig. 1 is for the cases of the sensor viewing 
at the center (6 =1.02°) with the solar zenith angles varying from 10°-80° at step 
of 10°. For comparison, the SeaWiFS results are plotted in the same figure. Fig. 1 
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shows that the implemented SeaWiFS atmospheric correction algorithm works as 
well for the MOS spectral bands as for SeaWiFS. We therefore conclude that, 
with appropriate computation of the Rayleigh scattering contribution at the MOS 
spectral bands, the current SeaWiFS atmospheric correction, with the lookup ta- 
bles of p a (A) + p ra (A), can be applied to MOS. 

2.3 A simple MOS destriping procedure 

The MOS radiance image has along-track stripes due to variations in the 
relative response of the individual detectors on the MOS CCD array (total of 384 
CCD detectors). Therefore, we have developed a simple destriping algorithm 
and applied it to the MOS radiance imageries. The MOS destriping procedure 
can be outlined as follows. First, for each scan (along the detector array) and a 
given spectral band, fit the radiance to a least-square cubic polynomial along the 
scan (the detector array) and compute relative gain at each detector (pixel), i.e., 

g{i,j) = X a„ i n / L(i,j), for i = 1 - 384, (5) 

n=0 

where L(i,j ) is the MOS measured-radiance for the detector number i and the 
scan number j for a given scene. Next, for each detector (pixel) select the median 
gain over all scans in the scene to derive the nominal gain factor for that detector, 
i.e., g{i) = Median[g(i,j )]. Finally, the MOS radiance image can be re-computed 
with the destriping correction L'(i,j) = g(i)L(i,j), where L' and L are the de- 
striped and original radiance, respectively. This simple procedure usually works 
quite well. Figs. 2(a) and 2(b) provide an example of results from the destriping 
algorithm for a MOS image acquired on February 28, 1998 in the Mediterranean 
Sea. Fig. 2(a) is the MOS original radiance image (443 nm) in which the along- 
track stripes are clearly evident, while Fig. 2(b) shows the same image after the 
MOS destriping algorithm has been applied. The destriping algorithm works 
quite well in this case, removing most of the striping effects with no obvious loss 
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of image structure, i.e., the physical properties of the image are preserved. In 
general, the efficacy of the algorithm depends mainly on how well the radiances 
along the detector array can be fitted with the cubic polynomials. To apply the 
implemented SeaWiFS atmospheric correction algorithm to MOS, however, the 
destriping procedure is usually not necessary. Without destriping, the results 
from atmospheric correction will be somewhat degraded, as the radiance striping 
adds noise to the process. The MOS destriping procedure usually improves the 
retrieved ocean optical properties significantly. 

2.4 The results from MOS compared with SeaWiFS 

We applied the atmospheric correction to both MOS and SeaWiFS for co- 
located images and compared the retrieved ocean and atmospheric optical prop- 
erties. The MOS radiance image was first destriped to remove the detector 
variations within pixels. The implemented SeaWiFS atmospheric correction al- 
gorithm was then applied to the MOS imagery. Two MOS-SeaWiFS co-located 
images acquired on January 29, 1998 in the Atlantic ocean and February 28, 1998 
in the Mediterranean Sea were first tested. These two scenes, acquired one 
month apart, differ significantly in their ocean and atmospheric optical proper- 
ties. Fig. 2 shows the MOS radiance image (443 nm) which was acquired on Feb- 
ruary 28, 1998 at a location of about longitude 3° and latitude 38° in the Mediter- 
ranean Sea. Fig. 3 shows both the MOS destriped and SeaWiFS radiance images 
(443 nm) which were acquired on January 29, 1998 at a location of about longi- 
tude -32° and latitude 27° in the Atlantic ocean. In comparing the MOS retrieved 
ocean and atmospheric optical results with that of SeaWiFS, we found that (i) the 
MOS retrieved aerosol optical thickness at the NIR band was usually a factor of 
2-3 times higher than that of SeaWiFS; (ii) the MOS retrieved e(7,8) which char- 
acterizes the spectral variation of aerosol optical properties is unreasonably low; 
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and (iii) the MOS retrieved normalized water-leaving reflectances [p w (X)] N in the 
visible are significantly different from those of SeaWiFS. Tables 2(a) and 2(b) 
provide examples of comparison results for typical co-located MOS 10x10 (5x5 
for SeaWiFS) pixel regions retrieved from these two cases. The parameters in the 
tables were obtained by averaging over the retrieved single pixel values in the 
co-located area (MOS 10x10 and SeaWiFS 5x5). The selected parameters in Table 
2 are the normalized water-leaving reflectance [p w {X)] N for bands 1-4, the ratio of 

aerosol single scattering reflectance between band 7 and 8 e(7,8), and the re- 

0 

trieved aerosol optical thickness at band 8 x a (8). The differences in e(7,8) and t a (8) 
between MOS and SeaWiFS are not shown in Table 2. Direct comparison of these 
atmospheric quantities is not relevant, since they depend on the solar and view- 
ing geometry of the observation, and that geometry is different between the two 
sensors. Furthermore, since there is about 90 minutes difference between co- 
located MOS and SeaWiFS observations, the atmospheric conditions may have 
changed. However, the e(7,8) value should be ~1 for typical marine aerosols, and 
it should certainly be > 0.5. Obviously, the results from SeaWiFS are more rea- 
sonable. The retrieved [p M ,(A)] w from SeaWiFS indicates typical two different 
ocean optical properties from these two scenes. The scene of the Atlantic ocean 
represents a typical clear ocean region with chlorophyll concentration -0.1 
(mg/m 3 ), whereas the scene from the Mediterranean Sea is a kind of turbid ocean 
water with chlorophyll concentration -0.8-1.0 (mg/m 3 ). Since we are applying 
an identical atmospheric correction process to the two sets of measurements, the 
large discrepancy in the retrieved [p w ]„ values between the two sensors can 
probably be interpreted as a difference in sensor calibrations. It is therefore nec- 
essary to re-calibrate one sensor to the other, to allow for meaningful compari- 
sons of the retrieved ocean optical properties. 
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3. A Vicarious Intercalibration for MOS 

As discussed in section 2, the sensor-measured radiance at the TOA is de- 
scribed by Eq. (1). Essentially, the first four terms in Eq. (1) are contributions 
from the atmosphere and ocean surface, and the last term is contribution from 
the ocean. There are mainly two unknowns in Eq. (1) for ocean color remote 
sensing: the aerosol optical properties and the water-leaving reflectance in the 
visible. Therefore, if one has knowledge of the atmospheric aerosol and ocean 
optical properties, one can essentially predict the sensor-measured radiance at 
the TOA for the MOS wavelengths [16]. These computed radiances can then be 
used to vicariously re-calibrate the MOS bands. Due to differences in the orbits 
of MOS and SeaWiFS, measurements of the same geographic location will be 
about one and one half hours apart. Since the atmospheric conditions are likely 
to change over that period, we can not expect that the atmospheric properties 
measured by SeaWiFS are valid for the MOS observations. We therefore assume 
that the gain of the MOS 868 nm band is unchanged, thereby using the aerosol 
concentration from the MOS measurements, and only accept that the aerosol 
model determined by SeaWiFS is still valid. Next, by using the SeaWiFS re- 
trieved aerosol models we can predict the atmospheric effects in the MOS im- 
agery, i.e., the first three terms in Eq. (1). The whitecap radiance contribution can 
be estimated in the same way as SeaWiFS [4]. Finally, using the SeaWiFS re- 
trieved normalize water-leaving reflectance, [p w {X)] N , the water-leaving radiance 
at the TOA in the MOS imagery can be computed according to Eq. (1), and the 
gain coefficients for the MOS bands can be derived. To reduce the variation of 
the derived gain coefficients with various scans, multiple scans within the MOS 
scene can be used to obtain coefficient data and derive a best fit for the MOS 384 


detectors. 
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In summary, the intercalibration procedure can be outlined as follows: 

1. find a MOS and SeaWiFS co-located scene in which both the SeaWiFS 
TOA radiances and the retrieved normalized water-leaving reflectances 
are relatively uniform; 

2. retrieve the aerosol models and [p w (X)] N values from the SeaWiFS meas- 
urements for the corresponding MOS pixels; 

3. theoretically predict the MOS measured-radiances at the TOA from the 
SeaWiFS data; 

4. obtain gain coefficients for the MOS bands 1-7 for all MOS 384 CCD de- 
tectors; and 

5. use the derived gain coefficients from multiple scans within the imagery 
and fit the gain coefficients data with the least-square cubic polynomials. 

We have applied the recalibration procedure outlined above to the two MOS 
scenes acquired on January 29 and February 28, 1998. The MOS scene has a 384 
detector-scan with image size of 384x384. We derived the MOS gain coefficients 
for the MOS 384 detectors at every 5th scan, thereby providing a total of ~75 gain 
coefficients for every detector of the MOS scene. Figs. 4(a)-4(g) provide the de- 
rived gain coefficients for the MOS bands 1-7 from scenes acquired on January 29 
and February 28, 1998, while Fig. 4(h) shows the MOS derived aerosol optical 
thicknesses at band 8. To clearly see the differences of the derived gain coeffi- 
cients from the two different MOS cases, we only plotted 50 representative data for 
each case in the figures (there are total of ~2.8xl0 4 data for each case). Appar- 
ently, the derived gain coefficients for the MOS bands 1-6 have very similar val- 
ues in the two different cases, indicating that they are nearly independent of 
temporal and spatial variations. The derived gain coefficients for band 7, how- 
ever, are different in the two cases. It appears that the MOS band 7 performance 
is related to the atmospheric optical conditions (see Fig. 4(h)) and its gain ad- 
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justment is in opposite to other bands (gain coefficient > 1). For the MOS bands 
1-6, we fitted the derived re-calibration gain coefficients from both cases to a 
least-square cubic polynomial (dotted lines in Figs. 4(a)-4(f)), while individual 
fits were derived for the MOS band 7 for two different cases. Clearly, the MOS 
band re-calibration adjustments are significant and they strongly depend on the 
MOS detector number. For example, the MOS band 1 has a re-calibration gain 
coefficient of -0.90 for detector 1, while it is -0.98 for detector 384. The MOS 
band 3 has the most changes (except band 7) with gain coefficient of -0.81 for 
detector 1 and -0.89 for detector 384. Table 3 provides the derived MOS re- 
calibration gain coefficients fitted with the least-square cubic polynomial as 

G(X, i)=ic n (X) i n , for i = 1 - 384, (6) 

n=0 

where i is the MOS detector number and c n ( X) is the fitting coefficient of the cu- 
bic polynomial for order number n. The gain fitting coefficients for bands 1-6 in 
Table 3 were derived with the least-square cubic fitting from the two MOS 
scenes, while the two set of band 7 gain coefficients were derived, respectively, 
from the MOS scene acquired on January 29, 1998 in the Atlantic ocean and on 
February 28, 1998 in the Mediterranean Sea. Therefore, for a given MOS band, 
only 4 re-calibration coefficients are needed for the 384 detectors. 

4. Results and Discussions 

We applied the derived MOS gain coefficients as in Table 3 to the MOS 
measured-radiance at the TOA, and retrieved ocean and atmospheric optical 
properties for comparison with results from the SeaWiFS measurements. Since 
the derived MOS band 7 re-calibration gain coefficients depend on the atmos- 
pheric optical properties, we have modified the atmospheric correction algorithm 
such that the correction can also be operated using the MOS bands 6 and 8. The 
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water-leaving reflectance at the MOS band 6, however, is usually not negligible. 
We have assumed a constant value of the MOS band 6 of 0.0959% which 

corresponds to a normalized water-leaving radiance of 0.045 (mW/cm 2 pm sr). 
This value was used in all results reported in this paper when using bands 6 and 
8 in the MOS atmospheric corrections. 

Figs. 5 and 6 provide the histogram (%) for the retrieved ocean parameters 
from MOS and SeaWiFS for the case of January 29 and February 28, 1998 for 
various situations. Figs. 5(a)-5(d) are, respectively, the retrieved normalize wa- 
ter-leaving reflectances (%) for bands 1-4 for the MOS data acquired on January 
29, 1998 in the Atlantic ocean, while Figs. 6(a)-6(d) are for the case of February 28, 
1998 in the Mediterranean Sea. For comparison, the retrieved parameters with- 
out the MOS re-calibrations are plotted in the same figures. There are four cases 
in each figure: (i) results from the SeaWiFS measurements with the bands 7 and 8 
used in the atmospheric corrections, (ii) results from the MOS re-calibrated radi- 
ances with the bands 7 and 8 used in the corrections, (iii) same as in (ii) except 
that the MOS bands 6 and 8 were used in the corrections, and (iv) results from 
the MOS original radiance data with the bands 7 and 8 used in the corrections. 
Table 4 shows the total number of retrievals (pixels) contributing to the histo- 
gram plots in Figs. 5 and 6 for cases of (i) to (iii). The MOS has ~6 times more re- 
trievals from each co-located scene than SeaWiFS because of its high spatial 
resolution. Note that using the MOS bands 6 and 8 in the atmospheric correc- 
tions yields slightly more retrievals than when using bands 7 and 8. Figs. 5 and 6 
show that the vicarious calibration improves the agreement significantly. Tables 
5(a) and 5(b) give quantitative comparisons of the retrieved parameters between 
the MOS (after re-calibrations) and SeaWiFS for the peak values in the histo- 
grams as shown in Figs. 5 and 6. For comparison of the retrieved atmospheric 
optical parameters, the peak values e(7,8) in the histograms are also listed in the 
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tables. Obviously, the differences in the retrieved normalized water-leaving re- 
flectance between MOS and SeaWiFS are reduced tremendously. Also, the MOS 
retrieved e(7,8) values are now reasonable and very similar to the values from 
SeaWiFS. Note that for cases in which the MOS bands 7 and 8 were used for the 
atmospheric corrections, two different calibration gain coefficients were applied 
for the MOS 750 nm band for cases of January 29 and February 28, 1998. How- 
ever, when the MOS bands 6 and 8 were used in the corrections, the MOS band 7 
reflectance data were simply not used, thereby allowing a consistent set of re- 
calibration gain coefficients for the MOS bands 1-6 and 8 to be applied for both 
cases. Both Figs. 5 and 6 and Table 5 show that, though the results of using MOS 
bands 7 and 8 in the atmospheric corrections yield slightly better agreement with 
SeaWiFS, good results can be obtained by using the MOS bands 6 and 8 with an 
assumed constant water-leaving reflectance value at band 6. Since the MOS band 
7 gain coefficient depends on the atmospheric conditions, which are highly vari- 
able in time and space, using the bands 6 and 8 for the MOS atmospheric correc- 
tion is more practical. 

To further test the efficacy of the vicarious re-calibration approach, we have 
applied the MOS re-calibration gain coefficients, which were derived from Janu- 
ary 29 and February 28, 1998 data, to a MOS image acquired on September 24, 
1997 at a location of about longitude 13° and latitude 45° in the Adriatic Sea, and 
compared the results to those obtained from a co-located SeaWiFS image. For 
this test, the destriping algorithm was not applied and the MOS bands 6 and 8 
were used in the atmospheric corrections. Figs. 7(a)-7(d) provide the histogram 
(%) of the retrieved water-leaving reflectances (%) for bands 1-4 from the MOS 
measurements for comparison with the SeaWiFS. The results from the MOS 
original calibrations are plotted in the same figures. The total retrievals contrib- 
uted in each plot in Fig. 7 are 6.76xl0 3 and 3.46xl0 4 for the SeaWiFS and MOS 



14 


(bands 6 & 8) case, respectively. It is truly remarkable that, with a time differ- 
ence of 4-5 months from the recalibration scenes and a different geographic loca- 
tion, the MOS derived water-leaving reflectances are still in a good agreement 
with SeaWiFS. It certainly improves the MOS retrieval results significantly from 
the original calibration. 

Since a two-band ratio of the retrieved normalized water-leaving reflectance 
[p w (A)] w in the visible is used to infer the ocean near-surface optical properties, 
e.g., the chlorophyll concentration can be related to either a ratio of band 2 to 5 
([P w ( 2 )]/^/[Pw( 5 )]a- ) or band 3 to 5 ( [p» v (3)]jv/[p M ,(5)] JV ) [9, 17], we have compared 
the MOS retrieved ratio values (after re-calibrations) with SeaWiFS for the Janu- 
ary 29, February 28, 1998, and September 24, 1997 cases. Figs. 8(a)-8(f) show 
histograms (%) in the retrieved ratios of the normalized water-leaving reflectance 
between bands 2 and 5 and bands 3 and 5 for various scenes. Figs. 8(a), 8(c), and 
8(e) are, respectively, the MOS and SeaWiFS retrieved normalized water-leaving 
reflectance ratio between bands 2 and 5 for case of January 29, February 28, 1998, 
and September 24, 1997, while Figs. 8(b), 8(d), and 8(f) are results of reflectance 
ratio between bands 3 and 5. In generating these figures, the MOS bands 6 and 8 
were used in the atmospheric corrections. Fig. 8 shows that, after re-calibration, 
the MOS derived ratio of retrieved normalized water-leaving reflectance agrees 
well with that of SeaWiFS. Therefore, MOS should be able to obtain similar chlo- 
rophyll concentration results as SeaWiFS. 

Finally, to compare results of the spatial distributions in the retrieved nor- 
malized water-leaving reflectance between MOS (after re-calibrations) and 
SeaWiFS, Figs. 9(a)-9(c) provide color images of the MOS retrieved normalized 
water-leaving reflectance (%) at 443 nm compared with the SeaWiFS measure- 
ments for a scene acquired on January 29, 1998 in the Atlantic ocean, February 28, 
1998 in the Mediterranean Sea, and September 24, 1997 in the Adriatic Sea. In 
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generating these images, the MOS bands 6 and 8 were used in the atmospheric 
corrections for retrieval of the MOS normalized water-leaving reflectances. Fig. 9 
shows that, after re-calibrations, MOS has a very similar results as from SeaWiFS. 

5. Conclusions 

We demonstrate that it is possible and efficient to vicariously intercalibrate 
two different ocean color sensors. In this study, the SeaWiFS retrieved normal- 
ized water-leaving reflectance and aerosol models were used as "truth" to re- 
calibrate the MOS spectral bands. After MOS band re-calibrations, the differ- 
ences of retrieved normalized water-leaving reflectances between MOS and 
SeaWiFS are much reduced. The MOS retrieved e(7,8) values are much more 
reasonable and very similar to the SeaWiFS measurements after re-calibration. 
Since the MOS band 7 re-calibration coefficients depend on the atmospheric con- 
ditions, we modified the atmospheric correction algorithm such that the MOS 
bands 6 and 8 can also be used for the corrections. Therefore, a consistent gain 
coefficients for the MOS bands 1-6 and 8 can be used for various MOS scenes 
obtained at different times and locations. We show the efficacy of the vicarious 
calibration approach by applying the method to a MOS scene acquired 4-5 
months prior to the data used in deriving the gain coefficients. The MOS results 
are in reasonable agreement with SeaWiFS. With this vicarious calibration ap- 
proach, the retrieved results from different sensors can now be meaningfully 
compared and possibly merged. With the same procedure one can also re- 
calibrate satellite sensors using in situ ocean and atmospheric optical property 
measurements. The proposed vicarious calibration scheme is applicable to other 
ocean color sensors, e.g., Japan's Ocean Color and Temperature Sensor (OCTS) 
and the French Polarization and Directionality of the Earth's Reflectances 
(POLDER). 
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Figure Captions 

Figure 1 . The errors in retrieved water-leaving reflectance at the MOS spectral 
bands 443 and 520 nm compared with the SeaWiFS spectral band configuration 
using the implemented SeaWiFS atmospheric correction algorithm for aerosol 
model of Maritime with RH = 80% and for various solar zenith angles. 

Figure 2. Results from the MOS simple destriping algorithm for a MOS image 
acquired on February 28, 1998 in the Mediterranean Sea; (a) the MOS original ra- 
diance image at 443 nm and (b) after applying the destriping algorithm. 

Figure 3. The TOA radiance at 443 nm measured by MOS and SeaWiFS for case 
of January 29, 1998 in the Atlantic ocean. 

Figure 4. The derived gain coefficients for the MOS bands from scenes acquired 
on January 29 and February 28, 1998 for (a)-(g) the MOS bands 1-7 and (h) the 
retrieved aerosol optical thickness at the MOS band 868 nm. 

Figure 5. The histogram (%) of the MOS retrieved normalized water-leaving re- 
flectances (%) with and without re-calibrations in comparison with the SeaWiFS 
measurements for case of January 29, 1998 (Atlantic ocean) for (a)-(d) as for the 
bands 1-4. For cases of the MOS after re-calibrations, results using both the MOS 
bands 7 & 8 and bands 6 & 8 for the atmospheric corrections are presented. 

Figure 6. Same as in Figs. 5(a)-5(d) except both MOS and SeaWiFS images were 
acquired on February 28, 1998 in the Mediterranean Sea. 

Figure 7. Same as in Figs. 5(a)-5(d) and 6(a)-6(d) except both MOS and SeaWiFS 
images were acquired on September 24, 1997 in the Adriatic Sea. For results of 
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the MOS after re-calibrations, the MOS bands 6 and 8 were used in the atmos- 
pheric corrections. 

Figure 8. The histogram (%) of the MOS retrieved ratio of the normalized water- 
leaving reflectances between band 2 to 5 and band 3 to 5 in comparison with that 
of SeaWiFS for cases (a) and (b) January 29, 1998 in the Atlantic ocean; (c) and (d) 
February 28, 1998 in the Mediterranean sea; and (e) and (f) September 24, 1997 in 
the Adriatic sea. 

Figure 9. The retrieved normalized water-leaving reflectance image (443 nm) 
from MOS compared with SeaWiFS for a scene acquired on (a) January 29, 1998 
in the Atlantic ocean; (b) February 28, 1998 in the Mediterranean sea; and (c) 
September 24, 1997 in the Adriatic sea. The MOS bands 6 and 8 were used in the 
atmospheric corrections. 
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TABLE 1. MOS and SeaWiFS band center wavelengths. 


Band # 

MOS 

X (nm) 

SeaWiFS 

X (nm) 

Difference 

AX (nm) 

1 

408 

412 

-4 

2 

443 

443 

0 

3 

485 

490 

-5 

4 

520 

510 

10 

5 

570 

555 

15 

6 

685 

670 

15 

7 

750 

765 

-15 

8 

868 

865 

3 


TABLE 2(a). MOS retrieved parameters compared with SeaWiFS for a co- 
located MOS 10x10 pixels for case of January 29, 1998. The [p w ] N is in %. 


Parameters 

MOS 

SeaWiFS 

Difference (%) 


4.453 

2.634 

69.1 

[P.( 2)1* 

4.860 

2.224 

118.5 

[P„(3)]„ 

4.093 

1.557 

162.9 

[p*( 4)1* 

2.029 

0.881 

130.3 

e(7, 8) 

0.198 

1.015 

— 

T.(8) 

0.091 

0.029 

— 


TABLE 2(b). MOS retrieved parameters compared with SeaWiFS for a co- 
located MOS 10x10 pixels for case of February 28, 1998. The [p w ] N is in %. 


Parameters 

MOS 

SeaWiFS 

Difference (%) 

[P.(D1# 

3.325 

0.903 

268.2 

[P„( Vh 

3.839 

0.945 

306.2 

[p.( 3)]* 

3.620 

1.090 

232.1 

[P.(4)] w 

2.488 

0.834 

198.3 


1.194 


e(7, 8) 
* a (8) 


0.487 

0.095 


0.029 
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TABLE 3. The derived MOS gain coefficients as G(A,0 = ]Te < (A)r, for i = 1-384. 


w-0 


MOS 

X (nm) 

c„0> 

C,(X) 

C 2 (X) 

0,0.) 

408 

0.9029 

3.5xl0' 4 

-9-OxlO- 7 

1.4xl0' 9 

443 

0.8453 

3.8xl0- 4 

-7.0x10-7 

6.5xlO- 10 

485 

0.8097 

3.8xl0- 4 

-5.3x10-7 

2.1xlO- 10 

520 

0.8693 

1.7x10-4 

-4.9xl0- 8 

2.8xlO* 10 

570 

0.8701 

1.8x10-4 

2.2x10-7 

-4.6xlO- 10 

685 

0.9287 

7.6x10-4 

-2.6X10- 6 

3.5xl0- 9 

750+ 

1.3208 

-3.5x10-4 

-2.9x10-6 

7.4x10-9 

750t 

1.2287 

5.1x10-4 

-5.3x10-6 

-3. 2x1 O' 10 

868 

1.0000 

0.0 

0.0 

0.0 


^ From case of Jan. 29, 1998. t From case of Feb. 28, 1998. 


TABLE 4. Total number of retrievals contributed in Figs. 5 and 6 for three 
different cases. 


Case 


Total # of Retrievals 


SeaWiFS 

MOS (7 & 8) 

MOS (6 & 8) 

Jan. 29, 1998 

2.27x104 

1.37x105 

1.40x10 s 

Feb. 28, 1998 

2.24x104 

1.38x10 s 

1.41x10 s 
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Table 5(a). The MOS retrieved parameters compared with SeaWiFS after 
MOS band re-calibrations for the case of January 29, 1998. The [p w ] N is in %. 


Parameter 
(Peak value) 

SeaWiFS 

(7,8) 

MOS 

(7,8) 

Diff (%) 
(7,8) 

MOS 

(6,8) 

Diff (%) 
(6,8) 

[P„ (Ok 

2.51 

2.41 

-4.0 

2.26 

-10.0 

[pj m N 

2.11 

2.16 

2.4 

1.96 

-7.1 

[P„( 3)] w 

1.51 

1.66 

9.9 

1.51 

0.0 

[p.( 4)] w 

0.86 

0.91 

5.8 

0.71 

-7.0 

£(7,8) 

0.983 

1.009 

— 

— 

— 


Table 5(b). The MOS retrieved parameters compared with SeaWiFS after 
MOS band re-calibrations for the case of February 28, 1998. The [p w ] w is in %. 


Parameter 

(Peak value) 

SeaWiFS 

(7,8) 

MOS 

(7,8) 

Diff (%) 
(7,8) 

MOS 

(6,8) 

Diff (%) 
(6,8) 

[P„ (1)]* 

0.76 

0.81 

6.6 

0.91 

19.7 

[P w (2)]„ 

0.86 

0.86 

0.0 

0.91 

5.8 

[P w (3)]„ 

1.06 

0.91 

-14.2 

1.01 

-4.7 

[p w m N 

0.76 

0.76 

0.0 

0.86 

13.2 

£(7,8) 

1.159 

1.159 

— 

— 

— 
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Figure 7. 
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